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An ideal compressible fluid is considered, with an equilib- 
rium density being a given function of coordinates due to pres- 
ence of some static external forces. The slow flows in such sys- 
tem, which do not disturb the density, are investigated with 
the help of the Hamiltonian formalism. The equations of mo- 
tion of the system are derived for an arbitrary given topology 
of the vorticity field. The general form of the Lagrangian for 
frozen-in vortex lines is established. The local induction ap- 
proximation for motion of slender vortex filaments in several 
inhomogeneous physical models is studied. 

PACS: 47.15.Ki, 47.32. Cc, 47.37.+q, 52.30.Cv 



I. INTRODUCTION 

Hydrodynamic-type systems of equations are exten- 
sively employed for macroscopic description of physical 
phenomena in ordinary and superfiuid liquids, in gases, 
in plasmas, in other substances. In solving hydrodynamic 
problems, it is admissible in many cases to neglect all 
dissipative processes and use the ideal fluid approxima- 
tion, at least as the first step. With this approximation, 
a dynamic model describing flow is conservative. The 
Hamiltonian formalism is a convenient tool to deal with 
such systems which makes possible to consider 

in a universal way all nonlinear processes. A big num- 
ber of works is devoted to application of the Hamiltonian 
method in hydrodynamics (see, for instance, the reviews 
H , 0| and references therein) . 

One of the most important questions, permitted for 
a universal consideration in the frame of canonical for- 
malism, is the question about integrals of motion of a 
dynamic system. Accordingly to the theorem of Noether 
, H , each conservation law of a system is closely con- 
nected with a symmetry of the corresponding Lagrangian 
with respect to some one-parametric group of transfor- 
mations of dynamical variables. It is well known that the 
conservation laws for the energy, momentum, and for the 
angular momentum follow from the fundamental prop- 
erties of the space and time, namely from homogeneity 
of the time and from homogeneity and isotropy of the 
space. Due to these properties, shifts and rotations of a 
system do not change its Lagrangian. The characteristic 
feature of the hydrodynamic-type systems is that they 
possess, besides the indicated usual integrals of motion, 
also an infinite number of specific integrals of motion re- 
lated to the freezing-in property of canonical vorticity Q- 
Jl2| ]. The reason for this is a basic physical property of 
fluids, relabeling symmetry. For instance, in isentropic 



flows the circulation of the canonical momentum along 
any frozen-in closed contour is conserved. In usual non- 
relativistic hydrodynamics, where the canonical momen- 
tum coincides with the velocity, the given statement is 
known as the theorem of Kelvin about conservation of 
the velocity circulation |l3| , Q . 

Existence of an infinite number of integrals of motion 
influences strongly dynamical and statistical properties of 
liquid systems. This is the reason why a clarification of 
structure of conservation laws is very important, as well 
as search for such new parameterizations for dynamical 
variables, which take into account the integrals of motion 
more completely. In many cases, even when a dissipation 
is present but its level is low, it is still correct to speak 
about integrals of the corresponding conservative prob- 
lem, because values of some of them are conserved with a 
high accuracy, especially on an initial stage of the evolu- 
tion, while the system has not proceeded to a state where 
a role of dissipation is significant due to large gradients. 
Besides this, conservation laws in physical systems, as 
a rule, are associated with definite geometrical objects. 
Usage of these associations promotes understanding and 
vivid imagination of everything that happens. In hydro- 
dynamic models, the frozen-in vortex lines are such geo- 
metrical objects, so the present work is devoted to study 
of motion of vortex lines in spatially inhomogeneous sys- 
tems. 

Hydrodynamic equations describe, in particular, an in- 
teraction between "soft" degrees of freedom of a system 
- frozen-in vortices, and "hard" degrees of freedom - 
acoustic modes. The presence of soft degrees of free- 
dom is explained by the fact that equilibrium states of 
the fluid are highly degenerated due to the relabeling 
symmetry. Thus, no potential energy corresponds to soft 
degrees of freedom, unlike the hard degrees of freedom. 
Due to dominating effect of elastic potential energy, hard 
degrees of freedom behave, typically, like a set of weakly 
nonlinear oscillators. On the contrary, dynamics of soft 
degrees of freedom is not dominated by a potential en- 
ergy, and usually it is highly nonlinear. In a limit of slow 
flows, when a typical velocity of vortex structure motion 
is small in comparison with the sound speed, a dynamic 
regime is possible, in which the hard degrees of freedom, 
corresponding to deviations of fluid density p(r, t) from 
an equilibrium configuration po(v), are excited weakly. 
Then, completely neglecting the sound, in the homoge- 
neous case po — const one arrives at the models of in- 
compressible fluid. For dynamics of vortices in incom- 
pressible perfect fluid the so called formalism of vortex 
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lines has been developed recently ||- which takes 
into account the conservation of topology of the vortic- 
ity field p5[ , Application of this formalism allows 

one to deal with a partially integrated system, where the 
topology is fixed by the Cauchy invariant p3[ . In pro- 
posed description the frozen-in solcnoidal vorticity field 
is considered as a continuous distribution of the elemen- 
tary objects - vortex lines. Such formulation of inviscid 
hydrodynamics as the problem of vortex line motion has 
been very suitable for study of localized vortex structures 
like vortex filaments. Also, it seems to be an adequate 
approach to the problem of finite time singularity forma- 
tion in solutions of hydrodynamic equations Jl7| . 

The goal of the present work is to extend the vortex 
line formalism to the case when equilibrium density po(r) 
is a fixed nontrivial function of spatial coordinates due 
to a static influence of some external forces. Such sit- 
uation takes place in many physically important mod- 
els. For examples, it can be the gravitational force for 
a large mass of an isentropic gas, both in usual and in 
relativistic hydrodynamics, or it can be the condition of 
electrical neutrality for the electron fluid on a given back- 
ground of ion distribution in the model of electron mag- 
netohydrodynamics (EMHD). The theory developed can 
be also applied to tasks about long-scale dynamics of the 
quantized vortex filaments in a Bose- Einstein condensate 
placed into a trap of a sufficiently large size. The vortex 
line formalism seems to be a universal and adequate tool 
for investigation of slow inviscid flows in inhomogeneous 
systems. For instance, it makes possible, in a simple and 
standard way, to analyze qualitative behavior of vortices 
without detailed consideration of basic equations of mo- 
tion for the fluid. Therefore the proposed approach can 
have advantages over other methods when complicated 
systems will be studied. 

As a concrete result, the local induction approximation 
(LIA) in vortex dynamics will be analyzed for several spa- 
tially inhomogeneous physical systems, namely for Eule- 
rian compressible hydrodynamics in an external field, for 
EMHD, and for vortices in trapped Bose-Einstein con- 
densates. A new equation of vortex filament motion will 
be derived, which takes into account the inhomogene- 
ity of these systems [the Eq.(^)]. As to relativistic hy- 
drodynamics in a static gravitational field, the proposed 
method gives a more complicated LIA equation than the 
Eq. ( p3|), as it has been shown recently by the present au- 
thor |18fl. [See also the most recent paper |19| about dy- 
namics of an ultrarelativistic fluid in the flat anisotropic 
cosmological models of expanding Universe, where the 
formalism of vortex lines has been applied to systems 
with Hamiltonian functionals depending explicitly on the 
time variable, and the effect of nonstationary anisotropy 
of the space on vortex dynamics has been studied.] 

This paper is organized as follows. A short review of 
Lagrangian formalism for fluid media is given in Sec. II. 
It provides a basis for development in Sec. Ill of the 
vortex line formalism for spatially inhomogeneous sys- 
tems. Then in Sec. IV, the method developed is applied 



to derive approximate equations of motion for slender 
nonstretched vortex filaments in three above mentioned 
physical models. 



II. LAGRANGIAN FORMALISM FOR A FLUID 

From viewpoint of the Lagrangian formalism, the 
freezing-in property of the canonical vorticity is due to 
the special symmetry of the basic equations of ideal hy- 
drodynamics (|]- g], pjj ], p2[ . As known, the entire 
Lagrangian description of a motion of some continuous 
medium can be given by the three-dimensional (3D) map- 
ping r ~ x(a, t), which indicates the space coordinates of 
each medium point labeled by a label a = (01,02,03), 
at an arbitrary moment in time t. The labeling a can 
be chosen in such a manner that the amount of matter 
in a small volume d 3 a in the label space is simply equal 
to this volume. With neglecting all dissipative processes, 
a dynamic model describing flow is conservative, so the 
equations of motion for the mapping x(a, t) follow from 
a variational principle 



SS = sJ £{x(a,<),x(a, t)}dt 



0, 



where the Lagrangian £ is a functional of x(a, t) 1 x(a, i), 
and also spatial derivatives. A very important circum- 
stance is related to the fluidity property of the media un- 
der consideration. The fluidity is manifested in the fact 
that the Lagrangian actually contains the dependence on 
x(a, t) and x(a, t) only through two Eulerian characteris- 
tics of the flow, namely through the field of density p(r, t) 
and the velocity field v(r,i), i.e., L = £{p, v}, with 



p(r, t) = det 



<9a(r,i) 



Or 



v(r,t) = x(a,i)| a=a(l%t) (1) 



Here a(r, t) is the inverse mapping with respect to x(a, t). 

A simple particular example is the Lagrangian of ordi- 
nary Eulerian isentropic hydrodynamics 



(2) 



C E uier = I ( Py - e{p) - pU (r) ) dr, 



where e(p) is the internal energy density, U(r) is the ex- 
ternal force potential, for instance, the gravitational po- 
tential. 

A less trivial example is the Lagrangian of relativis- 
tic isentropic hydrodynamics in a curved space-time 
with metric tensor gik(t,v) (i,k = 0..3; a, (3 — 1..3) 



2go a v a + g a pv a vP 



Here g = det || g^W is the determinant of the metric ten- 
sor, the expression in parenthesis is equal to the absolute 
value of the current four- vector n l — n{dx l /ds) pC) . A 
dependence £ (n) connects the relativistic density £ of 
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the fluid energy, measured in a locally co-moving refer- 
ence frame, with n. 

In plasma physics, the model of electron magnetohy- 
drodynamics is useful. EMHD follows in the limit of slow 
flows from the Lagrangian of electron fluid 



P- 



-p(A-v) 



(curl A) 

8^r~ 



dr, (3) 



where p(r, t) is the density of electron fluid, e is the elec- 
tric charge of electron, m is its mass, and c is the speed of 
light. The vector potential A(r, i) of the electromagnetic 
field determines the magnetic field B(r, t) by the relation 
B = curl A. In this paper, we will not need an explicit 
form of other terms indicated by the dots. 

The list of examples, of course, is not exhausted by 
three given models. All known hydrodynamic mod- 
els without dissipation, where the conservation of fluid 
amount takes place, can be described in this way. So the 
theory developed here is quite universal and applicable 
in various branches of physics where vortex phenomena 
occur. 

It follows from the definitions (Q) that dynamics of 
the density p(r, t) obeys the continuity equation in its 
standard form 



p t +V(pv) = 0. 



(4) 



The vanishing condition for variation of the action 
S = J £{p, v}dt, when the mapping x(a, t) is varied by 
5x(a,t), can be expressed in Eulerian representation as 
follows (the generalized Euler equation (TTJ ) 

<*-<-£)-*(£H(£)™ (5) 

This is merely the variational Euler-Lagrange equation 

d SC 5C 
dt (5x(a) (5x(a) 

for fluid particle dynamics. The equations (Q) and (||) 
determine completely evolution of hydrodynamic system. 

It was already mentioned that in all such sys- 
tems an infinite number of conservation laws ex- 
ists. The {p, v}-dependence means that the Lagrangian 
£{x(a, i), x(a, t)} admits the infinite-parametric symme- 
try group - it assumes the same value on any two map- 
pings xi (a, t) and X2 (a, t) , if they differ one from another 
only by some relabeling of the labels with unit Jacobian 

x 2 (a, t) = xi (a* (a) , t) , det 1 1 da*/da\ | = 1 . (6) 

Obviously, such mappings create the same density and 
velocity fields. According to the Noether's theorem |Q, 
[pj, every one-parametric sub-group of the relabeling 
group a* (a) with unit Jacobian corresponds to an inte- 
gral of motion. There are several classifications of these 
conservation laws. For instance, one can postulate that 
circulation of the canonical momentum p(r, t), 



P(r,t) 



5C 



5x(a, t) 



= i f scjp^y 

a=a(r,t) p \ <5v 



(7) 



along an arbitrary frozen-in closed contour 7(f) does not 
depend on time (the generalized theorem of Kelvin): 



f (P ■ dr) 

J-rit) 



const. 



We arrive at a different formulation, in terms of the so 
called Cauchy invariant, when consider the solenoidal 
field of the canonical vorticity fi(r, t), 



n(r, t) — curlp(r, t). 



(8) 



It is easy to check that application of the curl-operator 
to the equation (||) gives 



fl t = curl[v x n]. 
The formal solution of this equation is 



(9) 



n(r,t) = y<5(r-x(M))(n (a)Va)x(M)da, (10) 

where the solenoidal independent on time field l~io(a) is 
exactly the Cauchy invariant. The equation ( |l0| ) displays 
that lines of initial solenoidal field Jlo(a) are deformed 
in course of motion by the mapping x(a, tV keeping all 
the topological characteristics unchanged (T5J , |16[ . This 
feature of vortex lines dynamics is called the freezing-in 
property. 



III. HAMILTONIAN DYNAMICS OF VORTEX 
LINES 

To continue, it is more convenient to reformulate the 
problem in terms of density and canonical momentum. 
Let the system be specified by some Hamilton's func- 
tional H{p, p} 



n 



sc 



\S^ 



dr — C, 



(11) 



where the velocity v is expressed through the momentum 
p and through the density p with the help of Eq. ((?]). Let 
us note that the following equality takes place 



- I ( 5 J1 
p\5p 



(12) 



which is analogous to the formula (|7|). The Hamiltonian 
(non-canonical j^]) equations of motion for the fields of 
density and momentum follow from Eqs. (Q) and (|^). 
With taking into account the equality (|l^), they have 
the form (for detailed derivation see |ll|) 



Pt + V 



sn 



0. 



(13) 
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Pi 



5H 



curl p 

P 



SH 



(14) 



It is supposed in this work that the Hamiltonian has a 
minimum at some configuration {po(r), po(r)}. For sim- 
plicity, we will consider only systems without gyroscopic 
effects, i.e., po(r) = 0. Our purpose is to study slow 
flows near the equilibrium. In the regime under consid- 
eration, which corresponds formally to "prohibition" of 
excitation of the acoustic modes, the flow of fluid occurs 
in such a way that the density of each moving portion 
of fluid follows the given function po( r )- Therefore the 
equation (n3j) gives the condition 



(15) 



which means that after imposing the constrain p = p (r) 
the Hamiltonian does not depend anymore on the po- 
tential component of the canonical momentum field, it 
depends now only on the solenoidal component, i.e., ac- 
tually on the vorticity ft. It should be also noted that 
it is sufficient to take into consideration only quadratic 
part of the Hamiltonian, because the flow is supposed 
to be slow, so higher order terms, if any exists, may be 
neglected. Therefore, in further equations, H{tt} is actu- 
ally a quadratic functional of the vorticity field, though 
this fact will not be used in formal calculations. The 
condition (JL5|) implies validity of the formula 



SH 



curl 



SH 



(16) 



so the next equation for slow dynamics of the vorticity 
follows from ( ]14| ) 



curl 



curl 



SH\ ft 

sti) x ^) 



(17) 



This equation differs only by presence of the function 
Po{r) (instead of the unity) from the equation used in the 
works jlO| and Jl2j as a start point in the transition to 
the vortex line representation in homogeneous systems. 
Therefore all the further constructions will be done sim- 
ilarly to Ref. First, let us fix the topology of the 
vorticity field by means of the formula 

n(r,i) = J 5(t- R(a,t))(O (a)V a )R(a,i)da 

_ (fi (a)V a )R(a,i) 



a=R- 1 (r,t) 



(18) 



det||SR/0a|| 

where Jlo( a ) is the Cauchy invariant. The vector 

T(a,f) = (n (a)Va)R(a,t) (19) 

is directed along the vorticity field at the point r = 
R(a, t). It is necessary to stress that the information 



supplied by the mapping R(a, t) is not so full as the in- 
formation supplied by the purely Lagrangian mapping 
x(a, t). The role of the mapping R(a, t) is exhausted by 
a continuous deformation of the vortex lines of the initial 
field CIq. This means that the Jacobian 



J = dct||<9R/<9a|| 



(20) 



is not related directly to the density po(r), inasmuch as, 
unlike the mapping x(a,i), the new mapping R(a, t) is 
defined up to an arbitrary non-uniform shift along the 
vortex lines. Geometrical meaning of the representation 
( plf ) becomes more clear if instead of a we use a so called 
vortex line coordinate system (i/i(a), v% (a), £(a)), so that 
the 2D Lagrangian coordinate v = {u\,vx) £ AT is a label 
of vortex lines, which lies in some manifold TV, while a 
longitudinal coordinate £ parameterizes the vortex line. 
Locally, vortex line coordinate system exists for arbitrary 
topology of the vorticity field, but globally - only in the 
case when all the lines are closed. In the last case the 
equation ( [l8|) can be rewritten in the simple form 



n(r,t) 



' <j> S(r 



R(^e,i))R ? d£, (21) 



where R^ = <9R/<9£. The geometrical meaning of this 
formula is rather evident — the frozen-in vorticity field is 
presented as a continuous distribution of vortex lines. It 
is also clear that the choice of the longitudinal parameter 
is nonunique. This choice is determined exclusively by 
convenience for a particular task. Usage of the formula 

n t (r,t)=curl r J <S(r-R(a,t))[R t (a,i)xT(a,<)]da, (22) 

which follows immediately from Eq. (jl§), together with 
the general relationship between variational derivatives 
of an arbitrary functional F{£1} 



T x curl. 



SF{n{n}} 



(SR(a) 



(23) 



allow us to obtain the equation of motion for the mapping 
R(a, t) by substitution of the representation (jlj) into the 
equation (|T7j). As the result, dynamics of the mapping 
R(a, t) is determined by the equation 

[(n (a)V a )R(a) x R t (a)] p (K) = (24) 

It is not very difficult to check by a direct calculation 
that the given equation of motion for R(a, t) follows from 
the variational principle S J Cn dt = 0, where the La- 
grangian is 

£n = /'([R t xD(R)].(noVa)R)da-W{n{R}}, 

(25) 

with the vector function D(R) being related to the den- 
sity po( r ) by the equality 
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(Vr-d(R)) = p„(R) 



(26) 



For application to vortex filaments, the following form of 
the Lagrangian is more useful, where R = R(f, £, t): 



Cm = 



J d 2 v<£({R t x D(R)]-R ? )d£- H{fl{R}}. (27) 



It should be stressed that conservation in time of the 
fluid amount inside each closed frozen-in vortex surface 
is not imposed a priori as a constrain for the mapping 
R(a, t). All such quantities are conserved in the dynam- 
ical sense due to the symmetry of the Lagrangian ( p?| ) 
with respect to the group of relabelings of the labels v of 
vortex lines 



(28) 



Considering all one-parametrical subgroups of the given 
group of area-preserving transformations and applying 
the Noether's theorem m to the Lagrangian (|27]), it is 
possible to obtain the indicated integrals of motion in 
the next form (compare with Ref. j2lf| ) 

J* = J^{v u v 2 )d 2 i> Po (R)([Ri x R 2 ] -R^ (29) 

where ^(^i, v%) is an arbitrary function on the manifold 
M of labels, with the only condition $jayv = 0. 



IV. LOCAL INDUCTION APPROXIMATION 

When a particular task is being solved, the necessity 
always arises in making some simplifications. The vari- 
ational formulation for the dynamics of vortex lines al- 
lows us to introduce and control various approximations 
on the level of the Lagrangian (p7|), what in practice is 
more convenient and more simple than control of approx- 
imations made on the level of equations of motion. For 
example we will consider now the so called local induc- 
tion approximation (LIA) in dynamics of a slender non- 
stretched vortex filament. As known, in spatially homo- 
geneous systems the LIA yields an intcgrable equation 
which is gauge equivalent to the Nonlinear Schroedinger 
Equation [[22|. In general case, inhomogeneity destroys 
the integrability of LIA equation. Nevertheless, this does 
not reduce the value of LIA as a simplified model of fila- 
ment dynamics. 



A. LIA in Eulerian hydrodynamics 

At first, we will consider the Eulerian hydrodynamics, 
where the canonical momentum and the velocity coin- 
cide. Let the vorticity be concentrated in a quasi-one- 
dimensional structure, a vortex filament, with a typical 
longitudinal scale L being much larger than the width d 
of the filament. A typical scale of spatial inhomogeneity 



is supposed to be of order of L or larger. In such situa- 
tion, the kinetic energy of the fluid is concentrated in the 
vicinity of the filament, with the corresponding integral 
being logarithmically large on the parameter L/d. The 
LIA consist in the following simplifications. First, in the 
kinetic part of the Lagrangian (27), the dependence of 
the shape of vortex lines on the label v is neglected, i.e., 
the filament is considered as a single curve R(£, t). Af- 
ter integration over d 2 v the constant multiplier T appears 
now, which is the value of velocity circulation around the 
filament. Second, some significant simplifications may be 
done in the Hamiltonian. Generally speaking, exact ex- 
pression for the Hamiltonian implies derivation of the 
dependence v{po, ^1} from the following system of equa- 
tions: 

curl v = fl. div (po( r ) ■ v ) = 0, 

and subsequent substitution of v into the expression for 
the kinetic energy. After that one has to deal with a 
nonlocal Hamiltonian 



H 



{Po} 



' 1 ' G { a E po} {v^v^iv^p^dv^ 



Euler 2 

where the Green function G^' P °^( r ii r 2) has the follow- 
ing asymptotics at close arguments: 

Po(ri)<5 a /3 



Gi«' Po} (n,r 2 ) 



r 2 -> r x . 



afi 4 7 r|r 2 -r 1 | 
Therefore the Hamiltonian of a singular vortex filament, 



n 



d=0 



= i i a 



■<{E,Po} 
J af3 



(Ri.RaKAd&dfc, (30) 



where R' la — d^Ra^i) and so on, logarithmically di- 
verges. Taking into account the finite width d and the 
longitudinal scale L, it is possible to put, with a logarith- 
mic accuracy, the Hamiltonian of a thin vortex filament 
equal to the following expression 



where the constant A is 
.4 



— ln(- 



(31) 



(32) 



In accordance with the simplifications made above, the 
motion of a slender vortex filament in the spatially in- 
homogeneous system is described approximately by the 
equation 



[R 4 x R t ]p (R)/A = Vp (R) • |R C | -d 6 [ po(K) 



Rf 
IRf 



which is obtained by substitution of the Hamiltonian ( |3l| ) 
into the equation (E4T) . The given equation can be solved 
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with respect to R( and rewritten in terms of the geo- 
metrically invariant objects t, b, k, where t is the unit 
tangent vector on the curve, b is the unit binormal vec- 
tor, and k is the curvature of the line. As the result, we 
have the equation 



Kt/A = [V(lnpo(R)) x t] + «b, 



(33) 



applicability of which is not limited actually by the Eu- 
lerian hydrodynamics. Let us indicate at least two more 
physical models where the LIA equation (|33| ) is useful. 



quantized vortex filament |23] . At zero temperature this 
system is described approximately by the complex order 
parameter Q(r,t) (the wave function of the condensate), 
with the equation of motion (the Gross-Pitaevskii equa- 
tion) taking in dimcnsionless variables the form 



/<!>, = ( -|a + C/(r) - ^ + |*|" 



(34) 



where U(r) is an external potential, usually of the 
quadratic form 



B. LIA in EMHD 

The first model is EMHD, the Hamiltonian of which 
contains, besides the kinetic energy, also the energy of 
magnetic field B created by current of the electron fluid 
through the motionless inhomogeneous ion fluid. In prin- 
ciple, the Hamiltonian of EMHD is determined by the 
relations that follow from the Lagrangian C e , Eq. ([|) : 



curl v H B 

mc 



curl B 



4ne 



Po(r) • v, 



H 



EMHD 



f ( v 2 B 2 \ 



In spatially homogeneous system we would obtain the 
expression 



H 



EMHD 



P0_ 

8tt 



-g|ri-r 2 | 



ri -r 2 



-n(ri) • Jl(r 2 )dridr 2 , 



where the screening parameter q is determined by the 
relation 



q 2 = 



In inhomogeneous system q is a function of coordinates, 
with a typical value q. Let us suppose the inequalities 
qL ^> 1, and qd -C 1. One can see that the logarithmic 
integral analogous to the expression ([30j) is cut now not 
on the L, but on the skin depth A = X/q. Accordingly, 
for this case the constant A in LIA equation (^) is given 
by the expression 



r / mc 
Aemhd = -7— In 

471" 



We see that in ideal EMHD the LIA works better than 
in Eulerian hydrodynamics, due to the screening effect. 



U{t) = ax 2 + by 2 + cz 2 , 

and the constant p is the chemical potential. Let us sup- 
pose a > b > c. It is well known that the equation (|34| ) 
admits the hydrodynamical interpretation. The variables 
p and p arc defined by the relations 

p=|$| 2 , pp= ($V$-$V$)/2i. 

The corresponding Hamiltonian is 



K 



GP 



— J£ — 2 + (^( r ) - MP + y 



dr. 



In comparison with the ordinary Eulerian hydrodynam- 
ics, there is the term depending on the density gradient 
in this expression. However, with large values of the pa- 
rameter /i 2 /a, one may neglect that term in calculation 
of the equilibrium density inside the space region where 
the density is not exponentially small, and use the ap- 
proximate formula 



Po(r) 



U(r) 



if {Ai 2 >a; p-U{r)>0}. 



As known, the equation ([54|) admits solutions with 
quantized vortex filaments, the circulation around them 
being equal to 2tt. In these solutions, the density differs 
significantly from po( r ) only at close distances of order 
1/-JJ1 from the zero line. Far away, up to distances of 

order L ~ \J p/a S> 1/ <JJI, we have almost Eulerian flow. 
Therefore the LIA equation (|3^) is valid for description 
of slow motion of quantum vortex filament in trapped 
Bose-condensate of a relatively large size L, with the pa- 
rameter 



A = A 



GP 



(l/4)huy/a) 



The inequality p, 2 3> a ensures also the smallness of the 
filament velocity v/ ~ kAqp with respect to the speed 
of sound c s ~ y/JI, while the curvature of the filament is 
of order k ~ \fajp ■ 



C. LIA in Bose-Einstein condensate 

Another important physical model, where the equa- 
tion (33) may be applied, is the theory of Bose-Einstein 
condensate for a weakly nonideal trapped gas with a 



CONCLUSIONS 

Let us summarize briefly the main results of this pa- 
per. First, with neglecting acoustic degrees of freedom 



G 



in investigation of slow isentropic flows of a compress- 
ible perfect fluid in spatially inhomogeneous systems, 
the general form of variational principle for dynamics 
of frozen-in vortex lines has been found. The connec- 
tion from the basic Lagrangian given in terms of den- 
sity and velocity fields to the Hamiltonian of vortex lines 
has been provided, which allows one to analyze vortic- 
ity dynamics in complicated systems initially specified 
by the principle of least action. Second, this method 
has been applied to several physically important models, 
such as Eulcrian hydrodynamics in an external field, ideal 
electron magnetohydrodynamics on inhomogeneous ion 
background, and the Gross-Pitaevskii model for trapped 
Bose-Einstein condensate, in order to derive approximate 
equations of motion for vortex filaments. It has been es- 
tablished that a mathematical structure of the equations 
derived is the same in all the three cases, though param- 
eters have different physical meaning in each case. 

The final remark concerns possibility for development 
of an analogous approach in the general case, when acous- 
tic waves are important. Though at present moment this 
work has not been completed yet, but the author hopes 
it will be done in the future on the basis of equations 
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